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DIANA and selected applications* 

J. Fleischer, M. Tentyukov^ and O.V. Tarasov^ 

Universitat Bielefeld, Fakultat fiir Physik, Bielefeld, Germany 

New developments concerning the extension of the Feynman diagram analyzer DIANA are presented. We 
discuss new graphics facilities, different approaches to automation of momenta distribution and parallel processing 
facilities. Furthermore applications to ti production and Bhabha scattering are shortly discussed. 


The project called DIANA (Diagram ANA- 
lyzer)[l] for the evaluation of Feynman diagrams 
was started by our group some time ago. It was 
already used to calculate several processes [2]. 
The recent development * 1 of this project is do¬ 
cumented in this contribution. 

The pictorial representation of diagrams de¬ 
scribed in [3] includes three different kinds of 
postscript files. Now one more kind is available, 
the encapsulated postscript file containing parti¬ 
cle lines together with momenta flow [4] 2 : 



To use it, the user invokes the function 
\outInfoEPS(Fname,Height.Font.Fontsize). 
The picture will be saved in the file Fname. If 
Fontsize = 0 the particle images (e.g. e ± ) will 
not be printed. The width of the diagram will be 
calculated from “Height”. The diagram will be 
scaled to fit the EPS bounding box [0 0 Width 
Height]. 

Some parameters may be adjusted by the user. 
There are two functions, \setIEPSshift () and 
\setIEPSpars (), which may be invoked in the 
“initialization” section to reset default parame¬ 
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1 For details look at 

http: //www. physik. uni-bielefeld.de/'tentukov/diana.html 

2 For details see also http://www.physik.uni-bielefeld.de/ 

“tentukov/printing.html 


ters, see 

http: //www.physik.uni-bielefeld.de/ 

~tentukov/printing.html#OUTINFOEPS 

Particle identifiers can now be depicted by 
different fonts, sizes and colours. By default 
the particle image will be printed as in the 
corresponding entry of the propagator speci¬ 
fied in the ‘model file’. The user can set 
other images by means of the model extension, 
added in version 2.28, or by using the function 
\setparticleimage(A,image). The particle im¬ 
ages are produced by specific constructions, e.g.: 
W{y(10)+> —► W+. 

‘Blocks’ allow local shifting etc. of certain char¬ 
acters. A block is started by ‘{‘ followed by a 
keyword with parameters in parenthesis: 

{x(#) ... > - paints the content shifted 
along the abscissa (x). After the block, the cur¬ 
rent point is set to (old x, new y); 

-Cy (#) ... > - paints the content shifted 
along the ordinate (y). After the block, the cur¬ 
rent point is set to (new x, old y); 

{xy(#) (#) ... } - paints the content with 

shifts along both x and y. After the block, the 
current point is set to (old x, old y); 

{f(fontname)(#) ... } - sets font fontname 

scaled by =#= (in fractions of 1/10 of the current 
font size); 

{s(#) ... > - scales the current font by # 

(in fractions of 1/10 of the current font size); 

{c(#)(#)(#) ... > - sets RGB 3 colour. 

Each parameter must be a positive decimal num¬ 
ber in the region [0.. .1]. 

All sizes are in fractions of 1/10 of the base font 
size, i.e. “10” means the size of the current font, 

3 Red-Green-Blue - one of a standard colour model. 
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“20” is twice bigger, “5” is half size. 

As f ontname it is recommended to use only the 
standard 13 PostScript level 1 fonts: 

Courier Courier-Oblique 

Courier-Bold Courier-BoldOblique 

Helvetica Helvetica-Oblique 

HeIveti ca-BoId Helvetica-BoldOblique 

Times-Roman Times-Italic 

Times-Bold Times-Boldltalic 

£\|/|i|3oA, (Symbol) 

To set the image in the model, the user must 
add, after a line type, the character and the 
particle image. Example: 

[Wm,Wp;W;VP(num,ind:2,ind:1,vec,2);mmW; 

arrowWavy,3,2;W{y(10)-}] 

In this case the line will be of a photon-like ar¬ 
rowed curve of width 2 and amplitude 3, and the 
particle image will be W~. 

As alternative, the function 
\setparticleimage(particle,image) 
overwrites the image definition coming from the 
model. 

Topologies are represented in 
terms of ordered pairs of numbers like 
(fromvertex, tovertex) [1]. All external legs 
have negative numbers. 

Often the number of topologies is too large 
such that it is impossible to assign momenta 
to the lines in all topologies “by hand”. In 
that case, without special action, momenta are 
introduced automatically in an arbitrary man¬ 
ner (of course taking care of momenta conser¬ 
vation). The user names the loop momenta via 
the macro \loopmomenta in the “create” file, e.g. 
\loopmomenta(kl,k2,k3), and DIANA will as¬ 
sign momenta automatically using “AT”, “fc2” 
and “fc3” as the loop integration momenta. 

Sometimes it is important to keep some def¬ 
inite lines free from external momenta. If the 
users specifies SET _MARK_L00P=YES in the “cre¬ 
ate” file, the topology editor will be invoked in a 
special mode, and the user “clicks” (by mouse) 
which lines should carry bare integration mo¬ 
menta. All remaining momenta will be assigned 
automatically. 

Sometimes it is necessary to use more sophis¬ 
ticated distributions. For example, the user may 


want to use his favorite momenta like k — p \, k — 
P 2 ,k — P 3 , etc. assigned to some definite lines 
(see e.g. [5]). For such cases, DIANA provides 
the possibility to define momenta only for the 
virtual lines, and the full set of topologies (ob¬ 
tained by exchanging external momenta) will be 
defined from these “internal” topologies attach¬ 
ing the external legs. In this case momenta for 
internal parts are defined in terms of combina¬ 
tions of loop momenta (k in the above example) 
and some “abstract” tokens {pi,P 2 etc.), and for 
each topology DIANA expresses these tokens in 
terms of external momenta. 

Another example: occasionally topologies are 
generated by more complicated (“generic”) ones 
by scratching lines. In such cases one wants to 
stick to the momenta introduced for the lines 
which are kept [4] 4 . E.g. the user investigates 
the generic topology 
generictopology A = 

(-2,2)(-1,1)(1,3)(3,2)(2,4)(4,1)(3,4): 

pl+k2,pl+kl,kl,k2,k2-kl. 

Then DIANA will generate topologies from this 
one by scratching lines in the following manner: 



Topology A2_: 

3 



Topology A1_: 

_3 



4 For details look at http://www.physik.uni-bielefeld.de/ 

'tentukov/generictopology.html 
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The topology name is constructed by DIANA 
as follows: the name of the “generic” topology is 
appended by the index of “scratched” lines sepa¬ 
rated by underscore . 

To run an external command cmd, the function 
\exec(cmd) can be used. This function executes 
the command in the background without wait¬ 
ing for it to be completed. So, this function may 
be used to paralyse the evaluation of a process 
by running more then one FORM job simultane¬ 
ously. 

To avoid overloading the processor, by default 
at each time only one job is actually running while 
all the rest is waiting in the queue. On SMP 5 
computers, the number of simultaneously running 
jobs can by changed by means of a command line 
option -smp. Thus, -smp 8 tells DIANA to run 
on the computer 8 jobs simultaneously while all 
the rest is queued. 

To synchronize the TM 6 program with all 
started jobs, the function \waitall(timeout) is 
used. It suspends execution of the TM program 
until all jobs are completed. Each “timeout” mil¬ 
liseconds the function reports the number of jobs 
which are not finished yet. 

In case of a cluster of computers with a cur¬ 
rent directory shared by means of NFS 7 e.g., the 
function \exec() can use DIANA servers run¬ 
ning on other computers. To run DIANA as a 
server, the user must use the command line op¬ 
tion -d #, where # is the number of jobs which 
can run on this computer simultaneously. Thus, 
each computer on which the user has executed the 
command diana -d 1 -q is ready to perform the 
commands queued by the function \exec() (we 
assume that the current directory is shared by 
NFS among all computers in the cluster). Here 
-d 1 tells DIANA to run a daemon 8 accepting 
only one connection, and the option -q terminates 
the father DIANA process. For SMP computers 
the optimal argument for the -d option would be 
the number of processors. 

5 Symmetric Multiprocessing 

C TM is an abbreviation for “Text Manipulating”, see [1], 

'Network File System 

8 A program running in the background and listening to 

some port. 


Two functions, \exec() and \waitall(), per¬ 
mit the user to organize a simple parallel session 
for evaluations on SMP and/or cluster of inde¬ 
pendent computers with shared disk space. But, 
very often this is not enough. Indeed, let us sup¬ 
pose that all results must be collected into one 
resulting file log.all, while every job produces 
a file log.# where # is the order number of the 
job. Since jobs are running and completing in¬ 
dependently from each other, we can only collect 
all log.# into log.all after all the jobs are fin¬ 
ished. This leads to producing a lot of hies log. # 
at an intermediate step, which can overload a hie 
system. 

The simple solution is to allow some of newly 
started jobs to be synchronized with all previ¬ 
ously started jobs. Indeed, in this case after each 
job we can start another “slave” job, which ap¬ 
pends the log.# hie containing the result of the 
“master” job to the hie log.all. To obtain a 
correct order of the hie log .all, this copying job 
must be performed only after all earlier jobs are 
hnished. 

Another problem with cluster computations is 
that the optimal placement of the resulting hie 
log. # is usually a /tmp directory which is local 
against a current node 9 , on which a job “num¬ 
ber #” is performed. But to do this, the “slave” 
copying job has to “stick” to the “master” job, 
i.e. it must be performed on the same node as 
the ’’master” job. 

In order to create such a “slave” job, the func¬ 
tion \stick(cmd) must be used. This function 
is similar to \exec(cmd) with two exceptions. It 
performs the command cmd only after all earlier 
jobs are completed. This function in general is 
used to sum up all produced hies. That is why it 
performs cmd on the same computer as the pre¬ 
ceding job. 

As an example, let us consider the following 
script runf: 

#!diana -smp 1 -c runpar.tml 

\STARTSERVERS(phya25,phya26,phya27,phya28,phya29) 
\system(echo > log.all) 


9 DIANA assumes the “node” is the IP address of a server, 
so the conception of “nodes” is actually supported only for 
clusters. For SMP, the whole computer is assumed to be 
a single node. 
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\REPEAT(N) 

\exec(form -d i=\get(N) do.frm > /tmp/log.\get(N)) 
\stick(cat /tmp/log.\get(N) >> log.all) 

\stick(rm /tmp/log.\get(N)) 

\ENDREPEAT() 

\waitall(2000) 

The user enters: runf 186 200 , and the system 
executes: 

diana -smp 1 -c runpar.tml runf 186 200. 

The file runpar. tml contains definitions of var¬ 
ious TM functions and other settings; in partic¬ 
ular, it redefines the comment character as #. 
The macro \STARTSERVERS (list) checks if each 
server is working and, if not, it starts the server 
by means of the ssh. For example, for the host 
phya26 the following command will be performed: 

ssh phya26 cd CD ; diana -d 1 -q 

where “CD” is a current directory, e.g. 
/home/user/j obs. 

The operator \system(cmd) executes the com¬ 
mand cmd synchronously, i.e. it waits for the 
command to be completed and returns an exit 
code. Here it is used to produce an empty file 

“log.all”. 

All the instructions between 
\REPEAT(N) . . AENDREPEATO are cycled with 
N=186, . . . ,200. We assume that there is 
some folder file, say, tt. in with FORM in¬ 
put for each diagram. The FORM pro¬ 
gram do.frm evaluates a diagram by virtue 
of including a fold from the folder tt. in 
via an instruction like #include tt. in # n’i’. 
The macro definition i comes from the com¬ 
mand line form -d i=\get(N) do.frm, where 
\get(N) runs from 186 to 200. Each FORM 
job saves the result to the local directory, but 
the corresponding concatenation is performed by 
\stick(cat ...) on the same computer. At 
the end, all results will be collected in the file 
log.all, and all intermediate files \tmp\log.# 
will be removed. 

After all jobs are queued, the function 
\waitall(2000) will report every 2 seconds how 
many jobs are not yet completed. 

Recent applications of DIANA in our collab¬ 
oration are e + e _ annihilation into tt and Bhabha 


scattering. 

The tt production has been calculated includ¬ 
ing hard bremsstrahlung and various comparisons 
with other groups have been performed. For fur¬ 
ther information about the results we refer to re¬ 
cent notes on this subject [6,7]. 


Concerning Bhabha scattering, since we are 
heading for two-loop calculations, we perform all 
calculations in arbitrary dimension d = 4 — 2e. 
The algebraic part of this calculation has been 
presented in [8]. Here we discuss only the eva¬ 
luation of one-loop master integrals, in particular 
of box diagrams. The method applied uses eva¬ 
luation of difference equations [9]. We investigate 
now the following (scalar) box diagram with two 
(zero mass) photons in the s-cliannel: 


e 1 e" 

—►-•W/WWWWWV/WWVN.*-►— 


-*-•VWWWWVWWWWM 


The corresponding master integral in d dimen¬ 
sions can be represented as 
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where the ki(i = 1.. .6) are the following kine- 
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matical variables 


k\ = 1 — 


4?tt. 2 
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1 


k 2 = -4 mr(- + -- —-) = —m“z 

s t — 

s s 

^3 — 1 + 7 -;—, , k4 = 1 — 


fa = 1 


t — 4m 2 ’ 
t — 4m 2 


4m 2 


k 6 = 


4?n 2 


and $( 2 ,1, a) is the Lerch function 




y 

h k+a 


For s > 4to 2 and t < 0 it is easy to see that 
we remain in the analyticity domain of all occur¬ 
ring expressions, except for the factor (—s)fa~ 2 \ 
This result finally allows to expand the diagram 
in £ = 4=4. i n the two-loop calculation we need 
to expand up to order 0(e), which will yield fi¬ 
nite contributions when multiplied, e.g., by di¬ 
vergent counterterms, containing terms of order 
O(-). Due to the overall factor T(2 — |), coming 
from the infrared divergence of the diagram, the 
expansion of the various contributions in square 
brackets must be performed up to order 0(e 2 ). 
We also see that the imaginary part of the fi¬ 
nite part of the diagram comes from the above 
(—s)( 2 ~ 2 ) when expanded up to order O(e). 

Crossed diagrams can be obtained from this 
form as well if the result of the expansion yields 
analytic expressions (for which the analytic con¬ 
tinuation should be known in general) or again 
the kinematical variables are such that a nume¬ 
rical integration hits no singularities. 
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